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ABSTRACT 

The low-mass X-ray binary GS 1826—238 is presently unique for its consistently regular bursting 
behavior. In previous Rossi X-Ray Timing Explorer (RXTE) measurements between 1997 November 
and 2002 July, this source exhibited (nearly) limit-cycle bursts with recurrence times that decreased 
proportionately as the persistent flux increased. Here we report additional measurements of the burst 
recurrence time by RXTE, Chandra, and XMM-Newton, as well as observations of optical bursts. On 
a few occasions we measured burst recurrence times which deviated significantly from the earlier flux- 
recurrence time relationship, and most of these bursts occurred earlier than would be predicted based 
on the X-ray flux level. The epochs with early bursts were also accompanied by unusual broadband 
timing signatures, with the entire power spectrum shifting to higher frequencies. Concurrent XMM- 
Newton observations during one of these occasions, in 2003 April, indicate that an additional soft 
component may be present in the spectrum containing enough flux («30% of the total) to account 
for the burst recurrence time discrepancy. A self-consistent interpretation for the increase in soft flux 
and accompanying timing changes during 2003 April is that accretion disk extends down to smaller 
radial distances from the source than during the other observing epochs. The RXTE observations 
since 2003 April show that the spectral and timing properties have nearly returned to the previously 
established level. 

Subject headings: X-rays: binaries — X-rays: bursts — X-rays: individual (Ginga 1826-238) 



1. INTRODUCTION 

Type I X-ray bursts result from unstable thermonu- 
clear ignition of accreted material on the surfaces of 
weakly magneti c (B < 10 10 G) neutron stars (for a re- 
cent review, see IStrohmaver fc Bildstenl (|2006h ). Freshly 
accreted hydrogen and helium on the neutron star sur- 
face is hydrostatically compressed by new material at a 
rate m ~ 10 4 g cm -2 s _1 . In systems exhibiting bursts, 
the temperature and pressure at the base of the accreted 
layer slowly increase until the nuclear energy generation 
rate of the 3a-reaction becomes more sensitive to tem- 
perature perturbations than the radiative cooling. At 
this point the resulting thermonuclear instability leads 
to runaway burning of some or all of the matter that has 
been deposited since the previous burst. Typically, hours 
to days are required to accrete enough material to trigger 
the instability. 

There are close to 100 type I bursters that are known in 
the Galaxy 4 , and the vast majority are "atoll" type low- 
mass X-ray binaries (LMXBs) with luminosities above 
about 10 36 erg s _1 (|Hasinger fc van der Kllslll989f) . Al- 
though the basic physics of type I X-ray bursts is un- 
derstood, detailed comparisons between observations and 
theoretical models have had mixed success (e.g., Woosley 
et al. 2004). The most successful comparison has 
been with the "clocked burster" GS 1826-238 (also 
known as Ginga 1826-238), whose bursts are consistently 
quasi-periodic (|Ubertini et al.| [l999: Co cchi et al.l l200ll 
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ICornelisse et al.l [20031 ). Studies of the bursts from GS 
1826-238 began during the BeppoSAX mission, with the 
most detailed study of burst recurre n ce tim es and en- 
ergetics coming from iGallowav et al.l (|2004h . who ana- 
lyzed 24 bursts from the Rossi X-Ray Timing Explorer 
{RXTE) detected between 1997 November and 2002 
July. They found that the recurrence time decreased 
from 5.74 ± 0.13 hr to 3.56 ± 0.03 hr while the persis- 
tent (between burst) flux level increased by 66%, so that 
the recurrence time decreased almost precisely as 1/F X . 
Assuming F x cx M implies that the accumulated mass 
required for the instability to occur is approximately the 
same even as the accretion rate changes. The long burst 
durations (~100 s) and the low value of a (~40) - the 
ratio of the integrated bolometric persistent flux between 
bursts to the total bolometric burst fluence - both sug- 
gest that hydrogen c onstitutes a la rge portion of the 
fuel for these bursts (|Bildstenl [2000) . Further support 
for the H-rich fuel scenario comes from comparisons of 
the observed light curves with those predicted by time- 
dependent models, which also imply solar H and CNO 
composition in the material accreted by GS 1826—238 
(Heger et al. 2008). 

Additional RXTE observations have been made as 
part of a monitoring campaign from 2003 through 2007. 
These data reveal that the previously monotonic relation 
between burst recurrence time and persistent X-ray flux 
no longer fully describes the source behavior. Here we 
present an analysis of the full data set with the objective 
of more fully characterizing the complex burst behavior. 
We have also analyzed observations from the Chandra 
X-ray Observatory in 2002 July, and XMM-Newton in 
2003 April, both of which occurred simultaneously with 
RXTE, in order to also study the low-energy spectrum. 
In addition, we used optical observations made with the 
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UCT-CCD fast photometer at the South African Astro- 
nomical Observatory (SAAO) in 1998 June to determine 
the bursting behavior during gaps in the RXTE obser- 
vations. In § 2 we describe these observations and the 
details of our data analysis. We present our spectral and 
timing results in § 3. Finally, we discuss the implications 
of the results in § 4. 

2. OBSERVATIONS AND ANALYSIS 

2.1. RXTE 

In this paper, we utilized all RXTE observations of 
GS 1826—238, which has been observed at least once 
per year since 1997, except during 2001. Overall, there 
are more than 600 ks of good exposure time. We ana- 
lyzed data from the Proportional Counter Array (PCA; 
Jahoda et al. 2006) and the High Energy X-Ray Tim- 
ing Experiment (HEXTE; Rothschild et al. 1998) in- 
struments. The PCA has five identical co-aligned Pro- 
portional Counter Units (PCUs) sensitive to 2-60 keV 
photons. The HEXTE comprises two clusters, each con- 
taining four scintillation detectors sensitive to 15-250 
keV photons. Both instruments have large effective ar- 
eas (^6000 and 1400 cm 2 , respectively) and microsecond 
timing resolution. 

2.1.1. Recurrence Time Measurements 

Measurements of the burst recurrence times were made 
using PCA event mode data with 1 s bins, including all 
photon energies. The times of the bursts were defined to 
be when the flux exceeds 25% of the peak flux of each 
of the bursts. Sometimes bursts took place when the 
satellite was not taking any data (e.g., during passages 
through the South Atlantic Anomaly), and recurrence 
times were inferred if bursts occurred near enough to- 
gether to unambiguously know the number of intervening 
bursts. This method is acceptable because bursts from 
GS 1826—238 have never been observed to occur at irreg- 
ular intervals. If the time between two observed bursts is 
At, the recurrence time is estimated as At/(n+ 1), where 
n bursts are inferred in data gaps. Following Galloway et 
al. (2004), we adopted a fractional error of 2% on the in- 
dividual recurrence time measurements, and when bursts 
are inferred the error is scaled by l/(n+ 1) 1 / 2 because we 
assume n + 1 burst intervals in total. Although RXTE 
observations took place during 1999 and 2005, the dura- 
tion and spacing of these measurements did not allow for 
an unambiguous recurrence time measurement. Table 1 
shows the observation dates containing recurrence time 
measurements, the number of bursts that were observed, 
and the average burst recurrence time during each epoch. 
The observations taking place simultaneously (or nearly 
so) with RXTE are discussed below. 

2.1.2. Energy Spectra Analysis 

The evolving spectrum during the bursts was modeled 
using an absorbed blackbody, subtracting the pre-burst 
persistent emission as background. This approach is rel- 
atively stand ard for X-ray burst an alysis, and we refer 
the reader to iGallowav et al.l ((2004) for further details. 
In addition to the burst spectra, we studied the persistent 
emission by extracting color-color diagrams and photon 
energy spectra. In each case we excluded data from 500 
s before to 1500 s after a burst, the latter chosen to min- 
imize the residual blackbody flux seen during the decays 



TABLE 1 
Observation Log of GS 1826-238 



Date 






IxGC HIT . _1_ 11T1C 


(U.l.) 


Observatory 


h 

Bursts 


(hr) 


1997 Nov 5-6 


RXTE 


2 


5.88 c 


1998 Jun 7-12 


RXTE 


3 


5.67 ± 0.04 


1998 Jun 23-25 


SAAO/RXTE 


6/1 


5.13 ± 0.19 


2000 Jun 30-03 d 


RXTE 


8 


4.04 ± 0.06 


2000 Sep 24-27 


RXTE 


7 


4.11 ± 0.06 


2002 Jul 29-30 


Chandra/RXTE 


6/4 


3.54 ± 0.03 


2003 Apr 6-9 


XMM/RXTE 


16/7 


3.18 ± 0.14 


2004 Jul 19-21 


RXTE 


7 


3.47 ± 0.02 


2006 Aug 9-12 


RXTE 


11 


3.36 ± 0.07 


2007 Mar 8-10 


RXTE 


9 


3.53 ± 0.02 



Note. — Observations allowing unambiguous recurrence time mea- 
surements. Additional RXTE observations arc utilized in our study of 
the power spectra (sec § |2.1.3|) . 

a During three periods, simultaneous observations took place with 
RXTE. h All bursts observed by RXTE were also observed with the 
other observatory. Only one burst interval was measured during this 
epoch. This observation ended in July. 

of the bursts ([Thompson et~all[2005l) . Colors were cre- 
ated by extracting light curves for each PCA channel and 
observation, reading in the channel energy boundaries 
from the observation's response matrix, and interpolat- 
ing the counts spectrum to a standard grid. The soft 
and hard colors are defined as the counts ratio (3.5-6) 
keV/(2-3.5) keV and (9.7-16) keV/ (6-9.7) keV, respec- 
tively, and were obtained by integrating over the interpo- 
lated grid. To correct for the long-term drift of the gain 
of each PCA channel, the colors were normalized to the 
colors of the Crab calculated with the closest observation 
available to each GS 1826—238 observation. 

We extracted energy spectra from the persistent emis- 
sion between X-ray bursts using the standard software 
for RXTE data reduction (FTOOLS v.6.2). For the 
PCA, we used the "Standard 2" data, ignoring the first 
three channels and photon energies greater than 30 keV. 
HEXTE data below 25 keV were ignored, as were data 
above 100 keV due to poor statistics. The HEXTE chan- 
nels were rebinned so that there were a minimum of 
1000 counts per bin. To each observation (minus the 
times surrounding the bursts), we fit an absorbed dou- 
ble Comptonization model using XSPEC v. 11.3. This 
model was found by iThompson et alJ (|2005h to fit the 
broadband spectra of GS 1826—238 better than a single 
Comptonization model or empirical models like a cut-off 
power-law. For each measured RXTE burst recurrence 
time, the persistent flux was obtained from the best-fit 
model during the time since the previous burst. Al- 
though the spectral model only applies to photons from 
3-100 keV, we extrapolated the model from 0.1-3 keV 
and 100-200 keV to estimate the bolometric flux. 

2.1.3. Timing Analysis 

To relate the burst behavior to the broadband (per- 
sistent) timing properties, we also analyzed the rapid 
variability in the X-ray emission by producing a series 
of power density spectra (PSDs). For each RXTE obser- 
vation, including six observations that did not allow for 

5 Some of the fluxes had to be corrected for a 0?251 offset point- 
ing in R.A., which decreased the efficiency of the PCA and HEXTE 
by a factor of ~1.3. 
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a recurrence time measurement, we used 128 s segments 
of PCA event mode data with 2 -12 s bins (corresponding 
to a Nyquist freque ncy of 2048 Hz), and normalized after 
iLeahv et all {l983). As with the energy spectra and col- 
ors, we excluded data 500 s before to 1500 s after each 
burst. Separate groups of power spectra close in time 
were merged if no systematic differences in PSDs were 
observed. In this manner we obtained 18 power spectra. 
On average, each PSD contained 36 ks of accumulated 
data (with the smallest being 8 ks and the largest being 
83 ks). The contribution due to Poissonian statistics was 
estimated from the 1500-2048 Hz frequency band and re- 
moved (this also implicitly accounts for any decrease in 
power due to the PCA deadtime), and the spectra were 
converted to fractional rms squared. 

Qualitatively, each PSD can be characterized as having 
roughly equal power per decade in frequency, which ex- 
tends between a low and high frequency break. To quan- 
titatively describe each power spectrum, we fit a com- 
binati on of 4-6 Lorentzians following, e.g.. iBelloni et all 
( 2002), which overcomes the limitation of treating differ- 
ent power spectral components with intrinsically differ- 
ent models. Two to three of the Lorentzians are zero- 
centered and fit the broadband noise, giving a broad 
peak in the power times frequency iyP v ) representation 
between Vh at low frequencie s and vi and v u at high fre- 
quencies (jBelloni et al.l f2002) . There is also an additional 
component peaked at Vh to cover the "hump" typically 
seen near 2-3 Hz (but much higher during 2003 April, 
as shown below) in GS 1826—238. Additional narrow 
Lorentzians that are interpreted as quasi-periodic oscil- 
lations (QPOs; where the coherence factor Q = ^o/2A > 
2) were added if necessary to achieve acceptable fits. Due 
to limited statistics at high frequencies, only the lower 
frequency Lorentzian components (< 200 Hz) were use- 
ful for comparison. The frequencies are characterized 
using f ma x, the frequency at which the component con- 
tributes most of its variance per logarithmic frequency 
interval, which is equal to i/ maK = yJv^-V A 2 , where Vq 
is the centroid frequency and A is the half-width at half- 
maximum of the Lorentzian. Errors on the fit parameters 
were determined using A\ 2 = 2.71 corresponding to the 
90% confidence interval. 

2.2. Chandra 

The RXTE observations on 2002 July 29-30 oc- 
curred simultaneously with Chandra (also analyzed by 
iThompson et al.l [2005) , providing valuable low-energy 
coverage. A 68 ks observation was made with the Ad- 
vanced CCD Imaging Spectrometer (ACIS; Garmire et 
al. 2003), which is sensitive to photons from 0.3-10 
keV. High-resolution spectra were obtained by having the 
High Energy Transmission Grating (HETG; Canizares 
et al. 2005) placed in the optical path. The HETG 
spectrometer is composed of the Medium-Energy Grat- 
ing (MEG) and the High-Energy Grating (HEG). First- 
order MEG/HEG spectra were extracted with the Chan- 
dra Interactive Analysis of Observations (CIAO) v. 3. 4 
software. Responses matrices were generated using cali- 
bration v.3.4.0. Spectra were rebinned with at least 1000 
counts per spectral channel, giving statistically signifi- 
cant data for photon energies greater than 1 keV. 

2.3. XMM-Newton 



The RXTE observations during 2003 April 6-9 oc- 
curred simultaneously with an XMM-Newt on (Aschen- 
bach 2002) observation (also analyzed by iKong et al.l 
2007?) , lasting for 200 ks. Timing mode data was ac- 
quired from both the European Photon Imaging Camera 
(EPIC) and the Reflection Grating Spectrometer (RGS), 
although we only made use of data from the EPIC pn de- 
tector. The Optical Monitor was turned off for this ob- 
servation. Data analysis was performed with the XMM- 
Newton Science Analysis System (SAS) v. 7. 1.0. Photon 
energy spectra were extracted for energies between 0.5 
and 6 keV, and were rebinned with at least 1000 counts 
per spectral channel. Rather than using the "canned" 
response matrices, we generated one specific to our ob- 
servation. 

2.4. Optical 

During 1998 June, a set of three 7 ks RXTE observa- 
tions were made at ~1 day intervals, resulting in the de- 
tection of just one burst. At the same time, however, ob- 
servations of a small (50 x 33 arcsec 2 ) region surrounding 
the optical counterpart of GS 1826 —238 were made us - 
ing the UCT-CCD fast photometer (|Q'Donoghuel[l995T) . 
at the Cassegrain focus of the 1.9 m telescope at SAAO. 
On June 23 and 24, the observations lasted for 8 hours 
in an uninterrupted series of 5 s exposures; on June 25, 
the observation was 9 hours long using 2 s exposures. In 
each case, the coverage included (unambiguous) consec- 
utive optical bursts from GS 1826—238. We performed 
data reduction using the Image Reduction and Analysis 
Facility (IRAF), including photometry with the imple- 
mentation of DAOPHOT II (Stetson 1987). Point spread 
function (PSF) fitting was employed in order to obtain 
the best possible photometry, since there was moderate 
crowding of the counterpart and va riable seeing. The de - 
tails of this procedure are given bv lHomer efatl (|1998f ). 

Figure Q] shows the arrival times of the bursts seen in 
the optical and by RXTE. The bursts detected in the 
optical, coincident with (and immediately following) the 
burst detected by RXTE allowed an unambiguous mea- 
surement of the persistent flux and recurrence time dur- 
ing these observations. In addition, for the subsequent 
analysis, we associate the other two recurrence time mea- 
surements made possible by the pairs of bursts observed 
in the optical both before and after the burst detected 
by RXTE, with the X-ray flux measured during that ob- 
servation. 

3. RESULTS 

Within any individual observation, GS 1826—238 ex- 
hibited quasi-periodic bursting behavior, with the aver- 
age persistent flux in between bursts varying by at most 
a few percent. Figure [5] shows the variation of the burst 
recurrence time, burst fluence, and a as a function of the 
(absorbed) persistent X-ray flux using RXTE measure- 
ments. Except for the data from 2003 April, it is evi- 
dent that the recur rence time of the burs ts has a rough 
I /Fx dependence ([Galloway et al.1 12004!) . as has been 
observed previously, implying that typically, a similar 
amount of accreted fuel is required to trigger each burst. 
The typical scatter on the recurrence times within each 
epoch is just 1%. However, the 2003 April measurements 
are significant outliers, with recurrence times that are 
shorter than expected based on the persistent flux. Dur- 
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Fig. 1. — Burst arrival time plot during the optical observations 
in 1998 June. Intervening bursts (#3-5 & #8-10) were inferred to 
take place with consistent separations, although since the inferred 
times were during daylight at SAAO they could not be observed. 
The shaded bars show the observing windows for the nearly con- 
current RXTE observations. The dashed line assumes a constant 
burst recurrence time of 5.10 hr, although the burst behavior is 
more accurately described by a recurrence time of 4.92 hr for the 
first three observed bursts (1, 2, & 6, assuming constant intervals 
between 2 & 6), and 5.24 hr for the last four observed bursts (6, 7, 
11, & 12, assuming constant intervals between 7 & 11), using burst 
6 for each measurement. The apparent V-shaped residual plot is 
due to this change in recurrence time. If we assume the recurrence 
time changed discontinuously at burst 6, the residuals are only ~2 
min. 

ing this epoch, the recurrence time was about 3.2 hr even 
though the persistent flux was ^10% lower than in 2000 
June/September, when the recurrence time was about 
4.1 hr. Other measurements also deviate from the 1 /Fx 
trend, but to a lesser extent. For example, in 1998 June 
the recurrence time was typically near 5.65 hr, but about 
twelve days later optical burst intervals were 4.9 and 5.2 
hr even though the persistent X-ray flux (from RXTE 
observations) had not substantially changed. Therefore, 
the occasional discrepancy between the burst recurrence 
time and the expected value (given the persistent flux 
level and their previously observed relationship) seems 
to be a general property of GS 1826—238, and not some- 
thing that (for example) only occurs near the highest 
persistent flux levels. More recent RXTE data from 
2004, 2006, and 2007, have shown that the persistent flux 
versus recurrence time relationship has at other times 
obeyed the previously observed monotonic behavior. 

Figure [3] shows the fractional difference between the 
estimated bolometric flux and the flux derived from 
the empirical relationship between flux and burst recur- 
rence time (hereafter; referred to as the "recurrence time 
flux"), as a function of burst recurrence time. Note that 
this figure can be equivalently considered as the resid- 
uals of the x-axis in Fig. [2j top panel. Assuming the 
empirical curve describes the actual flux-recurrence time 
relationship, Fig. [3] suggests that the 2003 April X-ray 
flux measured by RXTE is underestimated by 30-40%. 

The discrepant measurements suggest either that the 
proportionality between the persistent 0.1-200 keV flux 
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Fig. 2. — Variation of the burst recurrence time, burst fluence, 
and a (= f Q t F p dt/ Jj^' F^dt) as a function of the absorbed per- 
sistent (between burst) flux in the range 0.1-200 keV (derived from 
extrapolating the model fit to RXTE data in the 3—100 keV range). 
Note that the unabsorbed fluxes are consistently ~5% larger. The 
dotted curve in the top panel shows the empirical fit to the data 
(excluding 2003 April), using At oc F~ 105 l IGallowav et alj|200lft ; 
the middle panel contains a linear fit. The 1998 June (opt) mea- 
surements in the top panel use optically-measured recurrence times 
and X-ray fluxes from the nearest RXTE observations (see Fig. [TJ. 
The lower panels only have one corresponding data point because 
RXTE measurements are only available for a single 1998 June (opt) 
burst. 

(derived from extrapolating the model fit in the 3-100 
keV range) and accretion rate was different than during 
other observations, or that the heating of the accreted 
layer (or some other factor) changes to alter the ignition 
depth at the same accretion rate. In order to discriminate 
between these two possibilities, we made a more detailed 
study of the burst characteristics, and of the persistent 
spectral and timing behavior of the source. 

3.1. The Energy Spectrum-Accretion Rate Relationship 

It is conventionally assumed that the accretion rate in 
LMXBs is proportional to the X-ray flux, and since we 
measure the flux in an instrumentally-determined pass- 
band, there is always the possibility that spectral varia- 
tions outside our passband will cause the proportionality 
to change between observation epochs. Here we explore 
this possibility in detail, both from the spectral shape 
(via the X-ray colors) within the RXTE passband (which 
is common to all the measurements), and in a broader 
band possible thanks to several occasions of contempo- 
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Fig. 3. — Fractional difference between the estimated bolometric 
flux and the "expected" flux derived from the empirical relationship 
between flux and burst recurrence time ( "recurrence time flux" ) as 
a function of the recurrence time. The symbols correspond to the 
same time periods as in Fig. [2] 




2.65 2.70 2.75 2.80 2.85 
Soft Color [(3.5-6) keV/(2-3.5) keV] 



Fig. 4.— Average RXTE/PCA colors of GS 1826-238 within 
0.2 days of the times of the bursts shown in Fig. [2] 

raneous Chandra and XMM-Newton observations. 

3.1.1. X-ray Colors 

Figure 2] shows the average RXTE/PCA colors within 
0.2 days of each burst. The individual 128 s colors (not 
shown here) show some scatter in the color-color diagram 
but all remain in the "island" state that typically char- 
acterizes bursting atoll LMXBs. Evidently, the colors 
from 1997, 2002, 2004, 2006, and 2007 are similar. The 
colors from 1998 and 2000 are comparable, although the 
soft color from 2000 and the hard color from 1998 are 
slightly smaller than the colors from the other periods. 
On the other hand, the colors from the 2003 April are 
both significantly smaller, with fractional changes in the 
soft and hard colors of about 4% and 3%, respectively. 

3.1.2. Broadband Spectra 

To test the extrapolation of our spectral model fits 
to energies below the PCA and HEXTE passband, we 
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Fig. 5. — Model spectra for the simultaneous Chandra/RXTE 
observations in 2002 July {dotted curve), and the simultaneous 
XMM-Newton/RXTE observations in 2003 April (solid curve). 
The Chandra/RXTE spectrum is restricted to E > 1 keV to reflect 
the energy range utilized in the fits; XMM-Newton provided good 
data down to 0.5 keV. 

also fitted the simultaneous 2002 July Chandra/RXTE 
and 2003 April XMM-Newton/RXTE broadband spec- 
tra 12.21 fc § 12. 3p . In both cases, two Comptonization 
components (comptt in XSPEC; Titarchuk 1994) were 
found necessary in ord er to accurately model the data 
([Thompson et al.ll2005l ). Each spectrum contains a com- 
ponent with ~0.1-0.2 keV seed photons, and a compo- 
nent with ~0.8-0.9 keV seed photons. One of these com- 
ponents is characterized by a ~6 keV electron plasma, 
and the other by a hotter electron plasma (~20 keV or 
greater), to fit the hard tail. We could not unambigu- 
ously associate the high or low seed photon temperatures 
with either the high or low plasma temperatures, how- 
ever, due to the model degeneracy that is introduced 
by the use of two Comptonization components. Table 2 
lists the best-fit parameter values for each spectrum and 
for the model coupling fcT s ,hot to fcTc.hot and kT s ^ co \d to 
fc?e,coid, plus fits to the 2003 April spectrum including 
an additional soft thermal component (see below). 

The model spectra for both epochs are presented in 
FigureEJ which shows that the major model-independent 
differences are that the 2002 July spectrum peaks (in 
vFy) at 1.3 keV, while the 2003 April spectrum peaks 
at 3.1 keV. RXTE observations without simultaneous 
low-energy coverage are, for the most part, insensitive 
to such changes in the soft flux, although the smaller 
soft colors (Fig. [4j during 2003 April are suggestive of 
this difference. The model fitting the broadband 2003 
April spectrum implies a 6% increase in the flux relative 
to the fit with only RXTE data, illustrating how a wider 
passband can reveal additional flux not apparent in the 
model extrapolated from the 3-100 keV band. Still, the 
discrepancy between the persistent flux and the "recur- 
rence time flux" is 30-40% (Fig. [3]). 

Relevant to the calculation of accurate bolometric 
fluxes for GS 1826-238 (see § [2X2]) is the question of 
additional soft spectral components. To test for the pres- 
ence of such components, we separately tried adding a 
multicolor disk blackbody (diskbb in XSPEC; Mitsuda 
et al. 1984) and a blackbody to our original spectral 
model, with the temperature of the soft component tied 
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TABLE 2 

Broadband Spectral Fits to the Persistent Emission 



Parameter 


2002 July 
2CTT a 


2CTT a 


2003 April 
2CTT+DBB b 


2CTT+BB C 


jV H (10 21 /cm 2 ) 
1. kT s (keV) 
1. kT c (keV) 

1. T 

2. kT s (keV) 
2. kT c (keV) 

2. T 

X? (d.o.f.) 


3.33 ± 0.01 
0.21 ± 0.01 
6.46 ± 0.03 
4.72 ± 0.01 

S5+ 4 - 65 
u - so -o.oi 

> 29.5 

1 fi5+ 893 
0.89 (1498) 


2.4 ± 0.1 
0.16 ± 0.01 

5.9 ± 0.4 

4.2 ± 0.1 
0.94 ± 0.02 

1 7 4+ 2 - 2 

n g+0.3 

1.22 (1152) 


3.73 ± 0.01 

[0.19 ± 0.01] 
5.9 ± 0.3 
4.4 ± 0.1 

0.75 ± 0.04 
I9.5i 2 ;* 
2.9 ± 0.2 

1.09 (1151) 


3.19 ± 0.01 
[0.15 ± 0.01] 
6.0 ± 0.3 
4.3 ± 0.1 
0.84 ± 0.05 

19.3^2^ 
2.8 ± 0.2 
1.02 (1151) 



Note. — Errors correspond to the 90% confidence interval for a single parameter. The 2002 July Chandra/ RXTE spectrum had statisti- 
cally significant data from 1—100 keV, and the 2003 April XMM-Newton/B.XTE spectrum had statistically significant data from 0.5—100 keV. A 
multiplicative constant was included in the models to account for differences in the absolute flux normalizations between the instruments. 
a XSPEC model: tbabs* (comptti + comptt2). Each comptt component assumes cylindrical geometry.* 3 XSPEC model: tbabs*(comptti + comptt2 
+ diskbb). The temperature of the disk blackbody at the inner disk radius was tied to the cooler seed photon temperature (indicated by "[ 
■ ■ ■ ]"). c XSPEC model: tbabs*(comptti + comptt2 + bbody). The blackbody temperature was tied to the cooler seed photon temperature (indicated 
by "I ■■■ ]»). 



to the seed photon temperature of the cooler Comp- 
tonization spectral component. We interpret the model 
with a multicolor disk blackbody as an approximation to 
the accretion disk spectrum, or partial covering of the 
accretion disk due to the disappearance of the Comp- 
tonizing medium over a limited volume, and the model 
with a blackbody as partial covering of the neutron star 
surface, or an extended optically thick boundary layer. 
In this picture, the soft photons acting as seed photons 
therefore emerge without being Comptonized. 

For the 2002 July spectrum, we found that the addi- 
tional component in each case (disk blackbody or black- 
body) did not improve the fit to the data, yet allowed 
for a 4-5% increase in the unabsorbed bolometric flux 
relative to the fit without the additional component. On 
the other hand, the fit to the 2003 April spectrum sig- 
nificantly improved in each case, with F-test probabil- 
ities of essentially zero that the improvement to the fit 
occurred by chance. The unfolded spectra and residu- 
als for the double Comptonization plus multicolor disk 
blackbody model is shown in Figure [HI (top & bottom 
panels). The residuals for the double Comptonization 
model derived using RXTE data alone, and using both 
XMM-Newton and RXTE data are also shown (middle 
panels). We stress, however, that the spectral fits for the 
models including a soft component should be considered 
illustrative. Accurately fitting a soft spectral component 
together with Comptonization components is problem- 
atic for larger seed photon temperatures because the 
comptt model assumes a Wien input spectrum rather 
than a thermal spectrum. This approximation, while 
leading to negligible changes to the spectrum at high 
energies, leads to an underestimation of the low-energy 
flux. 6 Compounding the issue is the interplay between 
the soft spectral component and the absorbing column 
density (see Fig. [5J inset); although the single-parameter 
uncertainty in the absorbing column density is very small 
(~0.2%), the two-parameter confidence contour between 
the absorbing column density and disk blackbody nor- 
malization (oci?i n ) suggests a much larger uncertainty. 

6 Nevertheless, the Comptonization model is more appropriate 
than empirical models such as a cutoff power law which, for T > 1, 
diverges at low energies. 



Because reliable interpretation of the fits is not possible, 
we do not include the soft component normalization val- 
ues in Table 2. Nevertheless, with the additional disk 
blackbody component the unabsorbed bolometric flux of 
the 2003 April spectrum is 3.7 xl0~ 9 erg cm~ 2 s _1 , 
which is ^50% larger than the values obtained with only 
RXTE. With an additional blackbody, the unabsorbed 
flux is 3.3 xl0~ 9 erg cm~ 2 s _1 , or ~30% larger than 
the values obtained with RXTE alone. Although the 
flux-recurrence time relation (Fig. [2]) used the absorbed 
fluxes, a direct comparison can be made by noting that 
the unabsorbed fluxes for the other observing epochs are 
consistently ~5% larger (because a constant absorbing 
column density was used when extrapolating the 3-100 
keV models to the wider 0.1-200 keV passband). The 
large fraction of soft flux contained in these components 
can clearly be seen in Figure [7l Therefore, despite the 
limited ability to acquire physical understanding of a soft 
thermal component in the 2003 April spectrum, it is a 
distinct possibility that the presence of one can account 
for the apparent disparity between measured persistent 
flux and the flux level that is expected given the burst 
recurrence time (Fig. 

3.2. Rapid Variability 

Each of the 18 average power spectra show very 
similar characteristic frequencies and rms amplitudes 
for the Lorentzian components representing the band- 
limited noise. The one major exception is the power 
spectrum from 2003 April, in which the characteris- 
tic frequencies are significantly higher. The 1998 June 
optical/ RXTE data show frequencies higher by a factor 
of ^2. A single low frequency QPO is present in seven 
of the PSDs, and in all but two cases this QPO appears 
to be associated with Vh . Five of the PSDs have two low 
frequency QPOs, and in each case the higher (lower) fre- 
quency QPO also appears to be associated with (v^). 
Similar low frequency QPOs have been observed in black 
holes and other neutron stars (e.g., Olive et al. 1998, 
Nowak 2000, Jonker et al. 200 0), and even in p revious 
observations of GS 1826-238 (|Barret et al.ll2000D . 

Three of the power spectra are displayed in Figure [SJ 
The top panel shows the power spectrum from 1998 June 
(optical) , the middle panel from 2003 April, and the bot- 
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GS 1826-238: 2003 April 




60 | 1 ' 

50 r ^ 

40 r \ + 2CTT (model derived from RXTE data only) 




4 : 2CTT+DBB: x', = 1 -09 |( , i j : 

^jWyiiii ^w^m 

-4 1 , , !_J L 

1 10 100 

Energy (keV) 

Fig. 6.— Unfolded 0.5-100 keV spectra and residuals for the 2003 
April observations using the double Comptonization plus multi- 
color disk blackbody model (2CTT+DBB, top and bottom panels). 
The two middle panels show the fit residuals for the double Comp- 
tonization model derived using RXTE data alone (upper middle 
panel), and using both XMM-Newton and RXTE data (lower mid- 
dle panel). The EPIC pn and HEXTE unfolded spectra are divided 
by 0.74 and 0.90, respectively, to account for the differences in the 
flux normalization between the instruments. The relatively large 
discrepancy between the P CA and EPIC norm alizations has been 
noted previously; see e.g., Yaqoob ct al. (2003). The residuals in 
the XMM-Newton data compared to the fit derived from RXTE 
data alone clearly indicate a substantial soft excess. By includ- 
ing XMM-Newton data, the double Comptonization model pro- 
vides an acceptable fit, but the fits with an additional soft thermal 
component are superior. The 2CTT+BB residuals are not shown, 
although they are nearly indistinguishable from the 2CTT+DBB 
residuals. For clarity, only every tenth residual for the XMM- 
Newton data are shown. Inset, top panel: Confidence contours 
(68% and 90% levels) between the absorbing column density and 
the inner disk radius of the multicolor disk blackbody component 
(assuming a 6 kpc source distance and 60° inclination). 

torn panel from 2006 August. The latter is more indica- 
tive of the timing characteristics seen during the other 
RXTE observation epochs. Figure [H shows the corre- 
lation between the two Lorentzian components that are 
present and well-constrained in all of the power spectra: 
Vb and i>h- The dotted line shows the best-fit linear trend. 
Similar correlations between the Lorentzian characteris- 
tic frequencies has been seen in many other sources (e.g., 
Psaltis et al. 1999, Wijnands & van der Klis 1999, Bel- 
loni et al. 2002, van Straaten et al. 2002), whether it be 
between the horizontal branch oscillation frequency and 
the lower frequency kHz QPO in Z-sources, or between 
the lower and upper frequency kHz QPOs in atolls, etc. 
Given that these correlated timing features are rather 
common in neutron stars and black holes, the present 
correlation does not come as a surprise. 




0.01 0.10 1.00 10.00 100.00 

Energy (keV) 



Fig. 7. — Unabsorbed 2003 April model spectrum containing an 
additional disk blackbody component or an additional blackbody 
component. The model is extrapolated to energies below the fitted 
range (0.5-100 keV). The dotted curves show the decomposition 
of the model components (the two Comptonized components are 
added together). The passband for RXTE is shown in the bottom 
of the figure, illustrating how a soft component could reside below 
RXTE coverage. Note how similar absorbed spectral shapes for 
E > 0.5 keV are achieved with a model containing a combination 
of either a blackbody and an absorbing column density of 3.2 X 
10 21 cm -2 , or a relatively brighter disk blackbody and a larger 
absorbing column density of 3.7 X 10 21 cm -2 . 

3.3. Early Ignition 

Although the X-ray spectral shape and timing charac- 
teristics are both atypical during 2003 April (and, to a 
lesser extent, 1998 June), it is also possible that some 
variation in the burst ignition properties contributes to 
the unexpectedly short bur st recurrence tim es measured 
during those observations. iBildsten ] ( j2000f ) pointed out 
that because the instability leading to bursts is a lo- 
cal phenomenon, the time for instability to develop de- 
pends on the accretion rate per unit area rather than 
the global accretion rate. A faster ignition time in 2003 
April could thus be explained by accretion over a more 
limited area of the surface of the neutron star relative to 
the other observing epochs. Even if the instability be- 
gins sooner, we expect the burst will consume all of the 
fuel that has been accreted since the previous burst be- 
cause, although the instability leading to the burst is a lo- 
cal phenomenon, the subsequent sprea ding burning front 
encompasses the entire stellar surface dSpitkovskv et al.l 
l2002t iBhattacharvva fc Strohmayerll2007f ). In this case, 
we would expect the burst fluences to be reduced approx- 
imately to the ratio of the recurrence times, i.e., ~3/4, 
assuming no change in the fuel composition at ignition, 
and a similar accretion rate in 2003 April as for the ob- 
servations at similar persistent flux levels. Instead, the 
fluences for the 2003 April are roughly similar to the 
bursts at higher or lower accretion rates (Fig. [2j middle 
panel) . 

One way the burst fluences in 2003 April could still 
match the measured values for the other bursts despite 
^25% less fuel is a change in the composition of the 
fuel. The nuclear energy generation rate during bursts 
is Qnuc = 1-6 + AX, where X is the mean hydrogen 
mass fraction in the layer, and we assume ~35% en- 
ergy loss due to neutrinos during the rp-process burn- 
ing (e.g., Fujimoto et al. 1987). A H-fraction that is 



1998 June (Opt) 




0.1 1.0 10.0 

Frequency (Hz) 



100.0 



Fig. 8. — Example power spectra from 1998 June (during the time 
of the optical bursts), 2003 April, and 2006 August, with the in- 
dividual Lorentzians superposed. The power spectra are presented 
in the vP v representation, which has the advantage of emphasizing 
higher frequency features and showing the maximum strength and 
frequency of the signal more directly. The components with char- 
acteristic frequencies of !4„ u^, v \ are labeled (the high frequency 
feature centered on i/ u is poorly constrained). The characteristic 
frequencies of the 2006 August PSD are more typical of the long- 
term average. The 2003 April PSD shows a dramatic increase in 
the characteristic frequencies, and the 1998 June PSD has moder- 
ately higher frequencies. 
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Fig. 9. — Characteristic frequencies of v b and v h for the 18 
average RXTE/PCA PSDs. The dotted line illustrates the correla- 
tion among these parameters. The additional symbols on the right 
(compared to Figs. 2-4) are from observations lacking a recurrence 
time measurement. 



TABLE 3 

Burst Profile Comparison 



Parameter 


2003 April 


2000 June/Sept 


>2002 July a 


T 


39.7 ± 0.6 


38.8 ± 1.2 


40 ± 4 


Tl 


18.0 ± 0.9 


17 ± 2 


19 ± 2 


T2 


47.3 ± 1.1 


43.5 ± 0.8 


45 ± 9 




24.5 ± 0.6 


26.5 ± 0.7 


25.5 ± 1.2 



Note. — The units of r, n, and T2 are seconds, and the units of 
-fpeak are 10 9 erg cm" 2 s — 1 . 
a Excluding 2003 April. 

«35% greater at ignition for the 2003 April bursts could 
thus compensate for the smaller accreted mass. Nat- 
urally, the mean H-fraction will be higher for bursts 
with shorter recurrence times, since less time has passed 
for steady H-burning to reduce the accreted composi- 
tion. The time to burn all the accreted hydrogen is 
ibum = 11(Z/0.02)- 1 (X /0.7) hr (e.g., Cumming & Bild- 
sten 2000), where Z is the CNO mass fraction, and X 
the accreted composition. Assuming the accreted hydro- 
gen begins burning immediately upon being accreted, the 
mean H-fraction in the fuel layer at ignition is then (for 
At < tbum): A (l - 0.5At/tbum). For X = 0.7 and 
Z = 0.02 this gives a H-fraction that is only ss6% larger 
for the 2003 April bursts. Therefore, changes to the com- 
position of the accreted material cannot explain how the 
fluences of the bursts in 2003 April are in agreement with 
the other observing epochs. 

Significant changes to the shape of the burst light 
curves may indicate unusual or unexpected nuclear burn- 
ing regimes. In order to characterize the burst profiles, 
we fitted a composite ex ponential curve t o eac h burst, 
and we refer the reader to iGallowav et alJ (|2006h for the 
details of the fitting procedure. We compared the 2003 
April burst decay parameters to those from the period of 
the most similar persistent flux (2000 June/September), 
and to the periods with the most similar burst recur- 
rence times (all observations after 2002 July, except 2003 
April). This exercise revealed that the shapes of the 
burst profiles in 2003 April are qualitatively similar to 
the other epochs (Galloway et al. 2004; Heger et al. 
2008). The parameter values are presented in Table 3. 
The only real indication of a significant difference be- 
tween the light curves is in the second decay constant 
(T2), which is somewhat longer for 2003 April than in 
2000 June/September. However, these parameters do 
systematically depend on the recurrence time due to the 
changing hydrogen fraction (see Fig. 4 from Heger et 
al. 2008). A more appropriate comparison might be 
with the bursts at the closest recurrence time, i.e., the 
bursts with a recurrence time of ~3.5 hr (2002 July and 
later, excluding 2003 April). Although this is a more in- 
homogeneous sample compared to 2003 April and 2000 
June/September, t-i for 2003 April is consistent with the 
others, though slightly higher (but not significantly) in 
the mean. There is a slightly longer second decay con- 
stant and lower peak flux for the 2003 April bursts com- 
pared to the bursts at the next nearest flux determined 
by RXTE, but, similar to the comparison of T2, it is con- 
sistent with the bursts at the nearest recurrence time. 

4. DISCUSSION 

During 2003 April, the recurrence times of thermonu- 
clear bursts from GS 1826—238 were unexpectedly short 
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given the persistent flux measured by RXTE and the 
previously-determined monotonic relation between these 
two parameters. At the same time, the spectrum from 
a simultaneous XMM-Newton observation (§ 13.1. 2[) sug- 
gests that an additional soft thermal component may be 
present at this time, which could increase the unabsorbed 
bolometric flux by up to 50%. A significant piece of ev- 
idence supporting a redistribution of the accretion en- 
ergy within the X-ray bands is the accompanying shift 
of the variability to faster timescales (Fig. 0, because 
correlations between higher power spectral frequencies 
and softer energy spectra are observed almost ubiqui- 
tously in black hole candidates and in many neutron 
stars. A direct relationship between QPO frequency and 
power-law in dex has been seen in many black hole can- 
didat es fe.g iVignarca et all 120031: iKajemci et all 12004 
120051 1200a iShaposhnikov fc Titarchukll2006D . and in at 
least five neutron stars (e.g., Titarchuk & Shaposhnikov 
2005). Moreover, episodes of rare soft thermal compo- 
nents (interpreted as the accretion disk) lasting for sev- 
eral months in the black hole candidate GRS 1758-258 
seem to be triggered by a decrease in the hard emission 
(jPottschmidt et aT1l2006h. Perha ps the most relevant ex- 
ample is from iFord et al.l (|1997t ). however, who showed 
that there is a direct correlation between the flux of the 
blackbody component and the QPO frequency in the 
burster 4U 0614+091. Although high frequency QPOs 
are not observed in the GS 1826—238 power spectra, the 
correlations between the frequencies of QPOs and broad- 
band noise components in other sources make it tempt- 
ing to speculate that the higher PSD frequencies in 2003 
April are indeed coupled with a soft spectral component, 
the presence of which is presumably related to changes 
in the accretion geometry. 

The correlations between faster variability and soft 
spectral components can be understood in the context 
of truncated disk models for accretion in X-ray binaries, 
which explain the low/hard to high/soft spectral transi- 
tions seen in both black hole and neutron star systems 
(for a recent review, see Done et al. 2008). For neutron 
star binaries at low L/L Edd (-0.04-0.08 in GS 1826-238, 
assuming d — 6 kpc), such models consist of a thin ac- 
cretion disk that is truncated at some radius far from the 
surface, probably due to evaporation of the accretion disk 
(jMever et al.ll200Ct iMaver fc Pringlell2007D . and an opti- 
cally thin 7 but geometrically thick hot inner flow that 
smoothly transitions to the boundary layer. Seed pho- 
tons from the surface of the star and from the inner edge 
of the disk cool the inner flow. At higher accretion rates, 
the inner radius of the disk moves inwards, which reduces 
the volume of the region occupied by the hot inner flow. 
Assuming that most of the variability comes from this re- 
gion, then a smaller volume would naturally account for 
the faster variability. In addition, the smaller disk radius 
will lead to a more prominent soft spectral component. 

One could argue that the soft photon flux should be 
present all the time, and perhaps just not detectable in 
the other observations due to poorer low-energy cover- 
age. According to truncated disk models, however, at 

7 Note that the spectral parameters in Table 2 do not support 
an optically thin interpretation. We defer a detailed analysis of 
the spectrum using more sophisticated models such as eqpair for 
a subsequent paper. 



low accretion rates the accretion disk is expected to be 
far from the surface of the neutron star (e.g., Esin et al. 
1997). If we assume that the power spectral break fre- 
quency v\, is associated with the inner disk radius and 
the orbital frequencies are Keplerian (i^Kep ex R~ 3 ^ 2 ), we 
expect the inner disk radius in 2002 July, for example, 
to be «3 times larger than in 2003 April, and the in- 
tegrated disk flux to be «3 times smaller (Frank et al. 
1992, p. 73). This is roughly consistent with the ratio of 
the absorbed soft fluxes shown in Fig. [5] 

Given that the 2003 April power spectral frequencies 
are the highest, and at the same time the burst recur- 
rence time suggests the accretion rate is the highest, it 
is worth considering that an alternative indicator for the 
mass accretion rate may be the timing data rather than 
the X-ray flux. For example, it has been proposed that 
t^b is positively correlated with the mass accretion rate 
(e.g., van der Klis 1994). If this were the case, however, 
we would expect a monotonic relationship between v\> 
and burst recurrence time (which we also assume to trace 
mass accretion rate). Such a relationship does not exist. 
To see this, one must only consider the data from 1998 
June (optical); the break frequency is the second high- 
est value even though the inferred accretion rate is the 
third lowest (only 1997 November and 1998 June have 
longer burst recurrence times, Fig. [2]). On the other 
hand, the two observing epochs showing the largest dis- 
parity between the recurrence time flux and the observed 
flux (Fig. ED, i.e., 1998 June (optical) and 2003 April, 
also exhibit the highest break frequencies. Therefore, it 
may be the case that 14, is a better tracer of the soft 
thermal flux and not the bolometric flux, although the 
absence of low-energy spectral coverage in the majority 
of GS 1826—238 observing epochs precludes a detailed 
comparison. 

An alternative explanation is that the bursts are ig- 
niting earlier, perhaps because the area over which the 
accretion occurs has decreased; however, any plausible 
mechanism to ignite the bursts earlier would also give rise 
to bursts that were significantly less energetic than oth- 
ers at comparable X-ray fluxes, whereas we find bursts 
in 2003 April with fluences comparable to those in other 
epochs. Thus, we conclude that the X-ray flux measured 
by RXTE in 2003 April underestimates the bolometric 
flux to a much greater degree than in other epochs, likely 
due to the presence of enhanced X-ray emission below —2 
keV. This implies that although the X-ray flux measured 
by RXTE in 2003 April was similar to what was observed 
in 2000 June/September, the accretion rate (and "true" 
bolometric flux) in 2003 April must have been at its 
highest level because the recurrence time was the small- 
est. This idea is consistent with truncated disk models, 
which explain the emergence of the soft spectral com- 
ponent by a reduction of the inner disk radius at high 
accretion rates. In this scenario, the previous and subse- 
quent RXTE observations of GS 1826-238 fortuitously 
maintained a fairly strict correspondence between X-ray 
flux in the RXTE passband and accretion rate due to a 
relatively constant spectral shape. Our results indicate 
that the spectral shape in GS 1826—238 starts to un- 
dergo a significant transformation when the bolometric 
flux is above about 3.5 x 10~ 9 erg cm -2 s _1 , or when the 
luminosity is above 1.5 x 10 37 erg s _1 (d/6 kpc) 2 . 
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There are other factors that could play a role during 
the periods of discrepant bursting behavior. One possi- 
bility is that the photon emission from GS 1826—238 may 
be anisotropic. Such emission is expected any time the 
scattering optical depth is much lower ove r a limited 
range of solid angles than in other directions ([King et al.l 
2001), and includes any source with an accretion disk. 
Recent work by Heger et al. (2008) suggests that emis- 
sion anisotropy does indeed influence the observational 
characteristics of GS 1826—238. In order for the burst 
timing and energetics from RXTE observations between 
1997 November to 2002 July to match theoretical igni- 
tion models, they found that £ p /£b = 1-55. Different 
anisotropy factors for the burst and persistent emission 
are expected because the radi ation can be att ributed to 
geometrically distinct regions (lFuiimotolll988h . If we as- 
sume that the burst emission is isotropic, the results of 



Heger et al. (2008) indicate that the persistent lumi- 
nosity may be underestimated by a factor of 1.55 in all 
observing epochs. To determine if anisotropic emission 
affected the 2003 April data to a higher degree than the 
other time periods, however, would require the results to 
be extended to include this epoch. 
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